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Summary. — People usually smile when astrophysicists assert that we are sons 
of the stars, but human life confirms this sentence: about 65% of the mass of our 
body is made up of oxygen, carbon occurs in all organic life and is the basis of 
organic chemistry, nitrogen is an essential part of amino acids and nucleic acids, 
calcium is a major component of our bones. Moreover, phosphorus plays a major 
role in biological molecules such as DNA and RNA (where the chemical codes of 
life is written) and our blood carries oxygen to tissues by means of the hemoglobin 
(an iron pigment of red blood cells). All these elements have been created in stars. 
I just list some examples related to human body, but also common element such 
as aluminum, nickel, gold, silver and lead come from a pristine generation of stars. 
The abundances in the Solar System are in fact due to the mixing of material 
ejected from stars that polluted the Universe in different epochs before the Sun 
formation, occurred about 5 billion years ago, after the gravitational contraction 
of the proto-solar cloud. Low mass AGB stars (1<M/M0<3) are among the most 
important polluters of the Milky Way, because of the strong winds eroding their 
chemically enriched envelopes. They are responsible for the nucleosynthesis of the 
main component of the cosmic s-elements. 

PACS 97.10.Cv - Stellar structure, interiors, evolution, nucleosynthesis, ages. 
PACS 97.10.Tk - Abundances, chemical composition. 



1. Introduction 

The majority of the isotopes heavier than iron (A>56) are synthesized by neutron 
capture processes. The observed heavy elements distribution shows the presence of two 
main components, correlated to different nucleosynthetic processes: the s (slow) process 
and the r (rapid) process (on the basis of the definitions given by [1] in their pioneering 
work). The r process requires high neutron densities, and it is believed to occur during 
explosive phases of stellar evolution (Novae, SuperNovae and/or X-rays binaries). The s 
process is characterized by a slow neutron capture with respect to the corresponding f3 
decay: stable isotopes capture neutrons, while the radioactive ones decay (/3 _ or (3 + ) or 
capture a free electron. These isotopes are mainly created in the Thermally Pulsing AGB 
(TP-AGB) phase of low mass stars (1 < M/M Q < 4), where freshly synthesized elements 
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Fig. 1. - Evolution of the positions in mass of the inner border of (top to bottom): convective 
envelope, H-burning shell and most energetic mesh of the He-burning shell, during the Thermally 
Pulsing AGB phase of a model with initial mass M=2M0 and solar metallicity.. 



are carried out to the surface by means of a recurrent mechanism called Third Dredge Up 
(TDU) (see [2] and references therein). In this phase the stellar structure consists of a 
partially degenerate carbon-oxygen core, an He shell separated from an H shell by the He- 
intershell region and by a convective envelope (see Fig.l). The energy required to supply 
the surface irradiation is mainly provided by the H burning shell, located just below the 
inner border of the convective envelope. This situation is recurrently interrupted by the 
growing up of thermonuclear runaways, driven by violent He-burning ignitions. As a 
consequence of a Thermal Pulse (TP), the region between the two shells (He-intershcll) 
becomes unstable against convection (for a short period), the external layers expand 
and, later on, the H shell burning temporarily dies down. In the He-intershell, He is 
partially converted into carbon. During the AGB phase, main neutron sources are the 
13 C(a,n) 16 reaction, active in radiative layers during the interpulse period [4], and the 
22 Ne(a,n) 25 Mg reaction, marginally activated within the convective shell originated by 
the TP. In order to obtain a sufficient amount of 13 C for the activation of the s-process, a 
diffusion of protons from the H-rich envelope into the 12 C-rich radiative zone is needed: 
the diffused proton are captured from the abundant carbon via the 12 C(p,7) 13 N(/3~) 13 C 
nuclear chain, leading to the formation of a tiny 13 C-pockct. 

2. The 13 C pocket and the nuclear network 

A major improvement in our stellar evolution code is the introduction of a physical 
algorithm for the treatment of the convective/radiative interface at the inner border of the 
convective envelope [3, 2]. During TDU episodes, the opacity of the envelope (H-rich) is 
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Fig. 2. - Abundances profiles in the 13 C pocket formed after the 3 th TDU episode (M = 2M Q 
model with solar metallicity). 



significantly larger than the opacity of the underlying H-exhausted (and He-rich) region. 
This fact causes an abrupt change of the temperature gradient at the inner border of 
the penetrating convective envelope. In this condition, the convective boundary becomes 
unstable, because any perturbation causing an excess of mixing immediately leads to an 
increase of the opacity and, in turn, to an increase of the temperature gradient. However, 
the steep pressure gradient that develops immediately below the formal border of the 
convective envelope limits the penetration of the instability, so that the average convective 
velocity should rapidly drop to zero. In order to mimic this behavior, we assume that 
in the region underlying the formal convective boundary, the average velocity follows an 
exponential decline, namely 




where d is the distance from the formal convective boundary, Vb ce is the velocity of the 
most internal convective mesh, Hp is the pressure scale height at the formal border of the 
convective envelope (defined by the Schwarzschild criterion) and (3 is a free parameter. 
The partial diffusion of protons in the top layers of the He shell gives naturally rise to 
the subsequent formation of a 13 C-rich tiny layer, the so-called 13 C pocket. In Fig. 2 we 
plot the chemical profiles of trace isotopes in the region around the 13 C pocket. The 13 C 
pocket (solid line) partially overlaps with a more external 14 N pocket (long-dashed line). 
The maximum of the 14 N coincides with the region where the protons diffused from the 
convective envelope at the epoch of the TDU were abundant enough to allow a full CN 
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cycle. This is also the region where the 23 Na (short-dashed-long-dashed line) has been 
efficiently produced by the 22 Ne(p,7) 23 Na reaction. Indeed, at the end of the previous 
TP, the mass fraction of 22 Ne (short-dashed line) in the top layer of the He-rich intershcll 
is of the order of 2xl0~ 2 and, around the 14 N maximum, this 22 Ne is almost completely 
converted into 23 Na, leading to the formation of a tiny 23 Na pocket. 
The introduction of the exponentially decaying profile of convective velocities automati- 
cally implies an assumption on the value of the free parameter f3: in the models presented 
in this work, we assume a value /3=0.1. In order to justify such a choice we perform 
some tests on a star with initial mass M = 2M Q and different metallicities (Z=Z Q and 
Z=1.0 x 10~ 4 ). We calculate the same sequences TP-Interpulse-TP, starting from the 
first TP followed by TDU, and we evaluate the consequences of changing this parameter 
in the range 0< (3 <0.2 (Cristallo et al., in preparation). While the amount of material 
interested by a single TDU episode increases by more than a factor of 3 from the case 
with no velocity profile ((3=0) and the extreme case ((3=0.2), the 13 C e //, available for 
the heavy elements nucleosynthesis, grows up with increasing (3, reaches a maximum 
in correspondence of (3=0.1 and then decreases down to negative values. We stress the 
fact that the net 13 C that contributes to the s-process nucleosynthesis is represented by 
13 C e //, defined as the difference between the 13 C and the 14 N mass fractions ( 14 N is 
in fact a strong neutron poison via the N(n,p) C reaction). Our result derives from 
a combination of two different physical processes: the growing of the 13 C e // is a direct 
consequence of the greater efficiency of the velocity profile algorithm, while its decreasing 
for large (3 values is due to the fact that the convection efficiency in the zone interested 
by the velocity profile starts to be too high. This leads the most external layers of the 
He-intershell to be fully mixed and to consequently show an envelope-like H abundance, 
instead of the H profile needed for the formation of the 13 C pocket. We therefore propose 
an efficient way to calibrate the free parameter affecting the velocity profile algorithm, 
confining its possible values in a narrow range, and we adjust it in order to obtain the 
maximum 13 C pocket expected for the s-process nucleosynthesis. 

The simultaneous solution of the stellar structure equations and a full network including 
all the relevant isotopes up to the termination point of the s-process path (Pb-Bi) requires 
a relevant computational power. For this reason, a post-process nucleosynthesis calcu- 
lation, based on AGB stellar models computed with a restricted nuclear network, was 
generally preferred [5]. The coupling of a stellar code with a complete nuclear network 
has not been feasible so far, but this limitation has been overcome thanks to last gen- 
eration of computers and to the adoption of smart algorithms to invert sparse matrixes 
of huge size. AGB models presented in this work are therefore calculated by using an 
extended set of nuclear processes including all chemical species involved in the s-process 
nucleosynthesis [6]. The models presented here have been obtained by including into 
the FRANEC stellar evolution code [7, 2] a full nuclear network, containing about 500 
isotopes (from H to the Pb-Bi-Po ending point) linked by more than 750 reactions. Re- 
action rates of isotopes involving charged particles are generally taken from the NACRE 
compilation [8], while the neutron captures cross sections arc mainly taken from [9]. 
Weak-interaction rates (electron captures, and (3 decays) are interpolated as a function 
of the temperature and electron density, while at temperatures lower than 10 6 K we 
assume a constant value equal to the terrestrial one. The presence of isomeric states of 
unthermalized isotopes, which lead to ramifications of the s-process flux, is followed in 
detail. This network is continuously upgraded according to the latest theoretical and 
experimental nuclear physics improvements. 
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Fig. 3. - Surface composition of a model with M=2M and solar metallicity after the fifth TDU 
(dotted curve) and the last TDU (solid curve). 



3. Low temperature C-enhanced opacities 

The FRANEC models benefit of another important improvement, whose effects are 
particularly important at low metallicities: the introduction of C-cnhanced low tempera- 
ture opacities. In principle, a stellar evolutionary code should account for the variations 
of the envelope chemical composition; in practice, only variations of the main constituents 
(H and He) are usually considered. Nevertheless, in the atmosphere of an AGB star, both 
C dredge-up and the possible conversion of C into N (by CN cycle) substantially affect 
the molecular contribution to the opacity for temperatures lower than 4000 K. This is due 
to the fact that, when a carbon excess is present (i.e. when the number ratio C/0>1), 
carbon bearing molecules (C2, CN, C2H2, C3) start forming and their contribute to 
opacity become dominant. Therefore, in order to overcome this limitation, we calculate 
low temperature C- and N-enhanced opacity tables, suitable for AGB model calculations 
[10]. The new AGB models succeed in reproducing the photometric and spectroscopic 
properties of their observational counterparts (see Section 6). 

4. Solar metallicity model 

The final C/O envelope ratio in the solar metallicity model is 1.88, while it becomes 
larger than 1 when the core mass is about 0.6 M Q and the bolometric magnitude is 
about —5.0 mag, in good agreement with the Galactic C-star luminosity function [11]. 
A 13 C pocket, whose mass extension decreases pulse after pulse, forms after each TDU 
episode. The 13 C in the first pocket is only partially burnt during the interpulse and 
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Fig. 4. - Surface compositions of models with M=2M and Z=6.0 x 10 3 (left Panel) and 
Z=3.0 x 1CT 3 (right Panel). 



the residual is engulfed into the convective zone generated by the subsequent TP. In this 
case, the 13 C burning takes place at the bottom of the convective shell (T ~ 1.6 x 10 8 
K) producing a very high neutron density, this fact implying interesting peculiarities on 
the following nucleosynthesis (see [12]). In particular, some branchings, which remain 
closed during a standard radiative 13 C burning and are marginally activated during the 
22 Ne(a,n) 25 Mg burning, are now open: in a very short temporal step (AT < 3 years) 
we obtain a consistent production of neutron-rich isotopes normally by-passed by the 
standard radiative 13 C s-process, among which 60 Fe, 86 Kr, 87 Rb and 96 Zr. Concerning 
the standard radiative s-process nucleosynthesis, all the elements (from Sr to Pb) result 
enhanced. We find that the abundance of Sr, Y and Zr at the first s peak, the so-called 
Is elements (light s elements), is comparable with the one of the hs elements (heavy s 
elements) Ba, La, Ce, Pr, Nd at the second s peak (see Fig. 3). Lead is under-produced 
with respect barium. The [hs/ls] (*) attained in the envelope when C/0=1 is in good 
agreement with those measured in galactic C(N) giants [13]. 

5. Intermediate metallicity models 

The 13 C(a,n) 16 reaction (the main source of neutrons) is primary-like, i.e. not 
directly affected by the metallicity of the pristine material. However, since the build-up 
of heavy nuclei requires neutron captures starting from Fe seeds, the s-process is expected 
to decline with metallicity (it is therefore of secondary nature). Nevertheless, while the 
iron seeds scale with the metallicity, when decreasing the number of iron seeds the number 
of neutrons available per seed is larger. The dependence of s-process yields at different 



C) We use the standard spectroscopic notation [A/B]=log(JV(A)/JV(B))-log(JV(A)/A/"(B)) Q . 
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Fig. 5. - Surface compositions of models with M=2M© and Z=1.0 x 10 3 (left Panel) and 
Z=1.0 x 1(T 4 (right Panel). 



mctallicities is therefore very complex and non linear (see e.g. [14]). At metallicities 
(3 -j- 6) x 10~ 3 the s-process distribution is peaked around the hs elements (see Fig. 
4). Moreover, we obtain an increase in the relative production of carbon and fluorine with 
decreasing the metallicity: this is a direct consequences of the greater efficiency of the 
TDU episodes, when the convective envelope penetrates into the 12 C-rich He-intershell 
(see [15]). We want to stress the fact that the production of 12 C is of primary nature (it 
is in fact produced by the activation of the triple a reaction in the He-intershell). The 
nature of the fluorine production, that basically depends on the nucleosynthesis of 15 N 
(see e.g. [16]) is more complex, being partly of primary nature (due to the contribution 
of the 13 C pocket) and partly of secondary nature (due to the contribution of the ashes 
of the H-burning shell). Finally, let us note that at these mctallicities the production of 
lead (mainly 208 Pb) starts being very efficient. 

6. Low metallicity models 

At low metallicities, most of the seeds are converted into 208 Pb, at the termination 
point of the s-process fluency (see e.g. [17]). Trends illustrated in Fig. 5 (which relates 
to the Z = lxlCP 3 and to the Z = lxlCP 4 models) confirm the previous sentence: a 
consistent lead production (Z=82) is in fact found in both models. As already stressed 
in the previous Section, the lower the metallicity is, the larger surface carbon overabun- 
dances are found: [C/Fe]~2.2 at Z=lxl(T 3 and [C/Fe]~3.0 at Z=lxl(T 4 . 
Low metallicity stars enriched in s-process elements show consistent enhancements in 
the hs region, spanning in the range 0.8<[hs/Fe]<2.3 . Even if our final [hs/Fe] value 
lays within the observed range, we cannot reproduce the observed spread with a single 
evolutionary model. In spite of this, we want to point out the importance of the opac- 
ity coefficients on the stellar evolution of low metallicity models. As already outlined 
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Fig. 6. - Surface compositions of a model with M=2M and Z=1.0 x 1CT 4 . The solid line refers 
to our standard model; the dashed line refers to a model with the same initial inputs but with 
a different opacity treatment (Z-fixed model). See text for details. 



in Section 3, in our models the opacity change caused by the variation of the internal 
chemical composition has been taken into account by linearly interpolating between ta- 
bles with different carbon, nitrogen and helium abudances [10]. However, the procedure 
followed by the majority of modelers is to adopt the opacity coefficients calculated at 
the initial mctallicity and to use them for the entire evolution. In order to evaluate the 
possible effects induced by these different opacity treatments, we compute another model 
at Z=lxl0~ 4 by using for the entire evolution the opacity coefficients calculated for the 
initial mctallicity (we refer to this calculation as the Z-fixed case) and we compare it 
with respect to the Z=lxl0~ 4 model already presented (standard case). The TP-AGB 
phase of the Z-fixed case is longer (a factor of 2 more than the standard model), this 
implying a natural increase in the number of TPs followed by TDU (N=49 with respect 
to the 15 found in the standard model). Consequently, the surface abundances of the 
Z-fixed case are larger with respect to the standard case (see Fig. 6), showing a better 
agreement with observations. We argue that these difference in the surface distributions 
are due to the mass loss rate, which results less efficient in the Z-fixed model with respect 
to the standard case, due to the fact that the surface temperature is always larger and 
the radius is less expanded [18]. As a consequence of so different mass-loss rate histories, 
in the Z-fixed case we find a large increase of the total dredged up material during the 
AGB phase with respect to the standard case and, therefore, the resulting final heavy 
element enhancements result larger. Our results seem therefore to suggest that the model 
computed with opacity tables at fixed mctallicity is more efficient in reproducing spec- 
troscopic data. There are however observational counterparts supporting the opposite 
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conclusion, driven by the study of C-stars in Dwarf Spheroidal Galaxies (DSGs); we refer 
in particular to the Draco Dwarf Spheroidal Galaxy [19, 20]. The observed giants be- 
longing to this DSph present surface temperatures in the range 3.55<logT<3.61 . Note 
that this range is never attained by the Z-fixed case, while it is entirely spanned by the 
standard case (see Fig. 4 in [10]). 

Another important aspect which has to be addressed is the choice of the mass-loss his- 
tory. Actually, observation cannot lead modelers in parameterizing the mass-loss at low 
metallicities, basically due to the fact that up to date a clear determination of the gas-to- 
dust ratio (needed to properly calculate the mass-loss) has not been feasible at non-solar 
metallicities [21, 22]. For this reason we adopt the same mass-loss history at all metallic- 
ities (see [2] for a description of our choice). Preliminary results show however that when 
using a different (milder) mass-loss [23] the surface overabundances increase on average 
for 0.2 deX (Cristallo et al. in preparation). 

Another source of uncertainty in modelling AGB stars at low metallicity is the choice of 
the mixing length parameter, which has been calibrated by calculating a Standard Solar 
Model. If using initial abundances from [24], we fit the observed properties of the Sun 
by using a value of the mixing length parameter a m .i =2.1 (see [25]). Due to the current 
lack of indications for the assumption of a different a m .i. value at low metallicities, we 
adopt the same value at all metallicities: preliminary results indicate however that when 
using a lower value (i.e. a TO .;.=1.8) the surface overabundances decrease on average for 
0.2 deX (Cristallo et al, in preparation). 

7. — Very low metallicity models 

In the case of very metal poor AGBs, we find that, for a given metallicity, there exists 
a lower mass limit for which normal AGB evolution may occur [10]. Stars with smaller 
initial mass experience protons engulfmcnt episodes during the first fully developed TP 
(see e.g. [26, 27]), which lead to a peculiar s-process nucleosynthesis, a low 12 C/ 13 C 
and a significant synthesis of primary N. Owing to the low entropy barrier, protons are 
engulfed within the He-intershell, where the high temperature and the large C abundance 
induce a violent H-burning flash. The consequent synthesis of 13 C allows the activation 
of the 13 C(a,n) 16 reactions and a very large neutron density is attained. Thus, the 
production of neutron-rich nuclei, normally by-passed by a standard s-process, can occur 
in these stars. Later on, when the intershell cools down, a deep TDU occurs. Since in the 
envelope the C/O becomes grater than 1, the large N abundance favors the formation of 
CN molecules, thus inducing a further increase of the radiative opacities. This new set 
of very low metallicity models, that represents the current frontier of AGB modelling, is 
still under analysis. 

8. Conclusions 

The present work demonstrates that, nowadays, the computational power allows 
the coupling between a stellar evolutionary code and a full nuclear network. A dif- 
ferent treatment of the internal border of the convective envelope with respect to a 
bare Schwarzschild criterion allows the formation of the so-called 13 C pocket. For the 
first time in the literature, it has been possible to directly compare theoretical models 
and observational data both from a physical (luminosities, surface temperatures) and a 
chemical (light elements and heavy elements abundances) point of view . Moreover, we 
furnish a uniform set of yields (from hydrogen to lead) at different metallicities. The 
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importance of the adopted mass loss rate and of molecules contribution to opacity has 
been pointed out. While the first problem is still under analysis, the second one has been 
solved by using opacity tables with enhanced carbon and nitrogen abundances. Finally, 
we have addressed the importance of very low metallicity AGB models, by describing 
first preliminary (and promising) results. 

* * * 

I sincerely thanks Prof. O. Stranicro and Prof. R. Gallino for their invaluable ex- 
pertise and their deep knowledge of the arguments treated in this paper. Without a 
continuous interaction and ideas sharing with them I would not acquire the needed sci- 
entific uprightness to make research in a proper way. 
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